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ABSTRACT 

The inner product provides a conceptually and algorithmically simple method for cal- 
culating the comoving distance between two cosmological objects given their redshifts, 
right ascension and declination, and arbitrary constant curvature. The key to this is 
that just as a distance between two points 'on' the surface of the ordinary 2-sphere 
S 2 is simply an arc-length (angle multiplied by radius) in ordinary Euclidean 3-space 
£ 3 , the distance between two points 'on' a 3-sphere S 3 (a 3-hyperboloid Ti 3 ) is simply 
an 'arc-length' in Euclidean 4-space £ 4 (Minkowski 4-space A4 4 ), i.e. an 'hyper-angle' 
multiplied by the curvature radius of the 3-sphere (3-hyperboloid). 

Key words: cosmology: observations — cosmology: theory 
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1 INTRODUCTION 

Liske (2000) recently showed a new variant for calculating 
the comoving distance X12 between two objects at cosmolog- 
ical distances, given their redshifts and angular separation. 
However, given two celestial positions (in right ascension and 
decimation) ai,5i,i = 1,2, before X12 can be calculated via 
equation (7) of Liske (2000), it is first necessary to calculate 
the angle #12 (i.e. a in the notation of Liske 2000) between 
these two positions. 

This angle can be calculated either by converting the 
angles to points in Euclidean 3-space, taking the inner prod- 
uct (dot product) of these two points treated as vectors, and 
inverting the cosine relation between the inner product and 
612, or via a spherical trigonometry formula involving sines 
and cosines. 

However, the calculation of X12 is very closely analo- 
gous to the calculation of #12, apart from the addition of a 
dimension, a change in signature of the metric in the case 
of hyperbolic space (an 'open' universe), and multiplication 
by the curvature radius. The generalisation from the arc- 
length R812 in three-dimensional space to an 'arc-length' X12 
in four- dimensional space is conceptually straight-forward if 
the inner product is used. 

In other words, use of the inner product in Cartesian 
coordinates provides a method of obtaining X12, given the 
density parameter f2 m , the cosmological constant Qa, the 
Hubble constant Ho, and the redshifts, right ascensions and 
declinations of two objects, z[,ai,5i,i = 1,2, via the def- 
initions and equations in eq. (33). No other equations are 
required. 



The use of the inner product to calculate 612 [denoted 
a in the notation of Liske (2000)] in three dimensions is 
described in Sect. ^, and its generalisation to the calculation 
of X12 (in four dimensions) is derived in Sect. 0. 



2 THE DISTANCE BETWEEN TWO POINTS 
IN S 2 AS AN ARC-LENGTH IN £ 3 

Given two celestial positions in spherical polar (e.g. equato- 
rial) coordinates (ai,Si),i = 1,2, these can be converted to 
Cartesian coordinates 

Xi — R cos Si cos on 
yi — R cos Si sin at 

Wi — RsinSi (1) 

on the celestial sphere of arbitrary radius R (e.g. R = 1). 

The standard Euclidean inner product on the two vec- 
tors cii = (xi, yi, Wi), i = 1,2, can then be expressed either 
as 

(ai, ai) = X1X2 + 1/12/2 + W1W2 (2) 
or as 

(a±, a±) = R 2 cos 9i2- (3) 

Equations (1), (^) and (^) imply the value of 612- Hence, 
the length X12 of a geodesic in S 2 between a\ and 0,2 is 

X12 =R 812 =R cos" 1 [(0,1,0,2) /R 2 ]. (4) 

In words, a distance in S 2 is simply an arc-length in 1Z 3 . 



© 0000 RAS 



2 B. F. Roukema 



As long as two vectors are represented in Cartesian coordi- 
nates, the inner product and thereby the angle and distance 
can easily be calculated. 



3 THE DISTANCE BETWEEN TWO POINTS 
IN S 3 {Ti 3 ) AS AN 'ARC-LENGTH' IN £ 4 (M 4 ) 

How can this be generalised to one more dimension and to 
the hyperbolic ('open') case? 

3.1 Friedmann-Lemaitre-Robertson- Walker 
coordinates 

Let us first write the standard Friedmann-Lemaitre- 
Robertson- Walker metric as 

ds 2 = -c 2 dt 2 + a 2 (t)[d X 2 + E 2 ( X )(d(5 2 + cos 2 8da 2 )] (5) 
where the dimensionless curvature is written 



fi K = Sim + — 1 
k = sign(fi K ) = 0,±1 



s(x) 



R smh(x/R) 
X 

R sin(x/R) 



k = 
k = 
k = 



-1 


+ 1 



[cf. equations (1), (2) of Liske (2000)], 

(c/H )(\n K \)-°- 5 k = ±i 



R 



undefined 







(6) 



(T) 



(8) 



and the 'proper' distance x [ ec L- (14.2.21), Weinberg (1972)] 
is 



X. 



da 



Ho -'i/(x+*') ay/Sl m /a - fi K + Qao 2 



(9) 



From here on, only spatial sections (hypersurfaces) at 
constant cosmological time are considered, i.e. dt = and 
a(t) = 1. 

Note that x nas a length dimension in this paper (e.g. 
h" 1 Mpc), whereas length units (e.g. h _1 Mpc) are (pre- 
sumably) included in the scale factor a(t) by Liske (2000), 
in order that x i s dimensionless. This requires some (arbi- 
trary) choice in the length scale of a(t) by Liske (2000) for 
the flat case (k — 0) , whereas here, there is no need to define 
R for the flat case. 



3.2 Cartesian coordinates, the metric and the 
inner product in 4-D 

What is the meaning of R for k = ±1? 

Let us first define four-dimensional Cartesian coordi- 
nates so that the two points in comoving 3-space become 
two points di = (xi, yi, Zi, u>i) in a four-dimensional space, 
via 

Xi = E(Xi) cos <5 i COS Oti 

yi = S(x» ) cos Si sin on 
Zi = E(xi) sin Si 

( R cosh(xi/fl) 
Wi = < 

1 R cos( Xi /R) 



k = 
k = 
k = 



-1 


+ 1 



(10) 



for i = 1,2, and similarly for any arbitrary point (z, a, S) [cf. 
eq. (12.4) of Peebles (1993)]. 

The metric of this four-dimensional space is 



k (dx 2 + dy 2 + dz 2 ) 
dx 2 + dy 2 + dz 2 



dw 2 



±1 
0. 



(11) 



Why is this a useful choice of metric? 

For k = +1, it is simply the obvious (by induction from 



£ ) metric defining Euclidean 4-space £ , i.e. 



da; 2 + dy 2 + dz 2 + dw 2 . 



(12) 



The meaning of R for the case k = +1 then follows. 
A spatial hypersurface at constant cosmological time is a 
3-sphere 5 3 . The 2-sphere S 2 is normally thought of as em- 
bedded in Euclidean 3-space £ 3 , since this is conceptually 
easy, although there is no mathematical necessity for the 
embedding. Similarly, S 3 can be thought of as embedded in 
£ 4 , in which case it has a centre located in £ 4 , but external 
to 5 3 and a radius of size R. 

The existence of the centre does not contradict the 
Copernican principle: all points in the physical space S 3 are 
equidistant from the centre [according to the metric of £ 4 , 
eq. (p. 21)], which itself is located in the mostly non-physical 
space £ . 

In the hyperbolic case, k = — 1, the intuitive meaning 
of R is less obvious. If the absolute value had not been used 
in eq. (||), then R would have had an imaginary value. Pee- 
bles [discussion near eq. (12.4) of Peebles 1993] points out 
that substituting R by iR results in the required equations 
and relations for hyperbolic space. How can one imagine a 
'negatively curved' sphere with imaginary radius iR? The 
relations 



cos id — cosh 
srniO — — isinh(9. 



(13) 



provide a clue. The fact that cos i9 is real but sin i6 is imag- 
inary suggests that the symmetry between the four coordi- 
nates needs to broken. 

This is why the metric [eq. (|ll|)] shows that symmetry 
between the four coordinates is indeed broken. The metric 



da" 



-(dx + dy + dz 



dw 



(14) 



and is of course just the metric defining Minkowski 4-space 
A4 4 , familiar from special relativity, in particular in the two- 
dimensional version, M 2 . 

Unlike the case of special relativity, the full space here 
does not have physical meaning. On the contrary, just as 
the physical space S is only a subset of the mostly non- 
physical £ 4 , there is a particular subset of A4 4 which does 
have physical meaning. This subset can be referred to as 
a 3-hyperboloid Ti. 3 , and defined in an equation uniting S 3 
and Ti. 3 by first defining the inner product 



, k (X1X2 + 2/12/2 + Z1Z2) + w±w 2 k = _., 
(01,02) = < . . , n (15) 

[ 21X2 + 2/i2/2 + ziz 2 « = 0. 

This clearly satisfies (in £ 4 and in the domain of M 4 re- 
stricted to 111 > \/ x 2 + y 2 + z 2 ) the properties required of 
an inner product that 

Vai, a 2 , as g £ 4 (resp..M 4 ), Vc £ 1Z 
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Figure 1. The 2-plane X-W for the case k = +1, after the 
isometries / and g have been used to shift the two points ai and 
d2 to a'j and a' 2 respectively, i.e. to (0, R) and (x' 2 ,w' 2 ) (resp.) 
in the X—W 2-plane. The circle S 1 in this plane is part of the 
3-sphere S 3 defined by {a, a) = R? . The inner product ^a^a^ 

evaluates to (a^a 1 ^ = Rw' 2 = R 2 cos8. The arc-length R9 (arc 
shown as extra-thick curve) is R8 = Hcos -1 (w' 2 / R). 

{cai,a-2) = 0(01,02) 

(0,1+0,3,0,2) = (ai, a 2 ) + (03, 02} 

(01,02) = (a 2 ,ai) 

(ai,ai) > 0. (16) 

The 3-sphere S 3 and the 3-hyperboloid Ti 3 of the co- 
moving space are then defined 

{a,a) = R\ (17) 

i.e. all points in the four-dimensional space of norm ||a|| = 
R, where the norm induced by the inner product is 

|a| = V< a, a >, (18) 

b ut where ther e is again a restriction to the domain w > 
■J 'x 2 + y 2 + z 2 in the case of A4 4 . 

Flat 3-space is the 3-plane w = 0. 

The restriction to the domain w > x' 2 + y 2 + z 2 for 
the k = — 1 case is required in order to satisfy the stan- 
dard definitions of inner product, norm and metric as always 
non-negative. In special relativistic terminology, where w is 
the time variable, these definitions are extended so that the 
invariant interval may be either time-like (positive by the 
convention here) or space-like (negative here). 

In the present case, an extension is also needed, because 
although all the points of the space Ti 3 lie in the 'upper 
cone' and represent 'time-like' vectors, the small difference 
vectors which are of interest in integrating along this surface 
are in the complementary domain w < *J x 2 + y 2 + z 2 . If 
relativistic terminology is adopted here (keeping in mind 
that the 'time' variable w here is purely artificial and has 
no physical meaning), then it can be said that although the 
invariant intervals (vector lengths) from the origin to the 




Figure 2. The 2-plane X-W for the case k = —1, after the 
isometries / and g have been used to shift the two points a\ 
and ct2 to a 1 , and a' 2 respectively, i.e. to (0, R) and (x' 2 ,w 2 ) 
(resp.) in the X-W 2-plane. The hyperbola H 1 in this plane is 
part of the 3-hyperboloid TH. 3 defined by (a, a) = R 2 . The in- 
ner product ^a'j^a^ evaluates to ^a'^a^ = Rw' 2 = i? 2 cosh0. 
The 'arc-length' R9 ('arc' shown as extra-thick curve) is R8 = 
Rcosir 1 (w' 2 /R). 

The 'hyper-angle' 6 is not an ordinary angle in this case. It can 
be thought as a parametrisation along the surface H 3 (or Ti 1 in 
the X—W plane) defined in such a way that R8 is a distance. 

points on Ti. 3 are all time-like, the difference element vectors 
tangent to the surface of Ti. 3 are all space-like. 

Hence, for k = ±1, let us replace equation (|ll]) by the 
more useful 'space-like' metric 

ds 2 |sp = k[k (dx 2 + dy 2 + dz 2 ) + dw 2 ] 

= dx 2 + dy 2 + dz 2 + k dw 2 (19) 

which can be integrated along Ti 3 (and, of course, also along 
S 3 ). 

3.3 The arc-length formula in four dimensions 

Does the definition in eq. (jfi|) lead to the following equation, 
which would appear to be the generalisation of eq. ([|)? 

{ficosh" 1 [(oi,o 2 ) /R 2 ] k = -1 

(ai - a 2 ,ai - a 2 ) k = (20) 

Rcos' 1 [(ai,a 2 ) /R 2 ] k = +1. 

3.3.1 Flat case 

The case k = is simply the Euclidean 3-distance in the 
3-plane w = 0. 

3.3.2 Curved cases 

The cases k = ±1 can be established by applying an appro- 
priate 'rotation' (isometry) and integrating the metric ds|sp 
between the 'rotated' positions a{ and a' 2 . 

As pointed out by Liske [just before eq. (6) of Liske 
2000], 'since the curvature of the three-dimensional space 
under consideration is constant, one can generate all totally 
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geodesic hypersurfaces from any given one by mere trans- 
lations and rotations', where the word 'rotations' is inter- 
preted loosely to include isometries of the 3-hyperboloid, 
H. 3 , as well as rotations of the 3-sphere, S 3 . 

So, either in the £ 4 representation of S 3 , or in the M 4 
representation of Ti. 3 , an isometry / can be chosen such that 
(0, 0, 0, 0) is kept as a fixed point and ai is shifted to 



f{ax) 
(0,0, 0,R). 



(21) 



After applying /, let us apply another isometry, a 3-rotation 
g in the x — y — z 3-plane, about (0, 0, 0, 0), which leaves w 
values unchanged and shifts a 2 to the point 



a' 2 = ff[/(a 2 )] 

= (x' 2 ,0,0,w' 2 ). 



(22) 



This leaves a{ unchanged, i.e. g(a[) — a[. 

Both a'i and a' 2 still lie on the surface defined by 
eq. (|17|), i.e. the relation 



kx 2 + w 2 = R 2 



(23) 



holds on this surface. 

To calculate the distance X12, let us integrate the metric 
ds|sp [eq. (19)] along the geodesic from a[ to a 2 (which lies 
in the X-W plane, see Figs [j], ^), i.e. 



X12 = 



d.s 



sp 



(24) 



and parametrise x and w in terms of a parameter 9 via 



cosh (w/R) k 
cos~ 1 (w/R) k 



-1 
+ 1 



(25) 



so that 



7?sinh6» 
R sin 9 



R cosh 9k- 
R cos 9 k - 

Ax 2 = w 2 d9 2 
dw 2 = x 2 d9 2 . 



-1 
+1 

-1 
+ 1 



(26) 



The endpoints of the integral a! x , a' 2 , then become 
h=0 (27) 



and 
9 2 = 



using eqs ((230, (E3). .-, 

The integral [eq. (|J)] is then 



cosh' 1 (w' 2 /R) k = -1 
cos' 1 (w' 2 /R) k = +1 



(28) 



X12 



ds 



sp 



+ k dw 2 



[using eq. (19)] 
r«3 



\/ w 2 + k x 2 d9 



[using eq. (26)] 



Rd9 



[using eq 
= R9 2 



Rcosh- 1 {w' 2 /R) k 
R cos" 1 (w' 2 /R) k 



[using eq. (g|)] 

i?cosh~ 1 [(a' 1 ,a , 2 ) / "R 2 
i?cos~ 1 [(ai,a 2 )/i? 2 ] 

[using eqs ©, (21), (22)] 



-1 
+ 1 



7?cosh~ 1 [(ai,a 2 ) / R 2 ] k = 



Rcos 1 [(ai, a 2 ) /R 2 ^ 
[since / and g axe isometries] 



k = 



-1 
+ 1 



-1 
+ 1 



(29) 



Thus, eq. (Efj) is the correct generalisation of eq. 



3.3.3 The flat case as a limit of the curved cases 

Why is the inner product used differently in the curved and 
flat cases in eq. (po[)? 

The reason can easily be seen by taking the limit as 
R — > 00, after isometries / and g have been applied as above. 

If k = ±1, but R max{i, y, z}, i.e. R 3> x' 2 , then 



X12 = 



7?cosh- 1 [(a' 1 ,a 2 ) /R 2 ] 
Rcos- 1 [(a' 1 ,a' 2 ) /R 2 ] 

i?cosh _1 (iu 2 /_R) k 
R cos~ 1 (w' 2 /R) k 

i?sinh _1 (x 2 /_R) k ■■ 
Rsm^ix'i/R) k : 

R (x' 2 /R) 
x' 2 , 



-1 
+ 1 



-1 

+1 

-1 

+1 



i.e. 



lim 

For k = 0, 



X12 



X 



12 = y/(x' 2 - 0) 2 + 2 + 2 = x' 2 = 



lim 



X12- 



(30) 
(31) 

(32) 

ii-*oo,fc = ±l 

Thus, the flat case is a limit of the curved cases as ex- 
pected. 

4 CONCLUSION 

So, just as a distance in S 2 is an arc- length in 1Z S , the dis- 
tance between two objects at cosmological distances in a 
curved universe can be thought of as the arc-length corre- 
sponding to an 'hyper-angle' in four- dimensional Euclidean 
or Minkowski space and is thus obtained directly from the 
inner product. 

This is algebraically equivalent to the solutions of Os- 
mer (1981), Peebles (1993), Peacock (1999) and Liske (2000) 
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but is expressed completely via the definitions and equations 
(6), ©, (1), (|), (10), © and ©. 

That is, the complete formulae for calculating an FLRW 
comoving distance between two objects at cosmological dis- 
tances, given Qra, ^a, Ho, {zi,Ui,8i),i = 1,2, are the fol- 
lowing: 

r i? sinh(x/.R) fe = -1 
E( X ) = { X k = 

[ i? sin(x/i?) fc = +1 



J? = 



X = 



(ai, 02) 



(c/iJo)(|n K |)-°- 5 fc = ±1 
undefined = 



c / da 

^° Ji/(i+z<) a^Q. m /a - Sl K + Q A a 2 
Q K = Q, m + f^A — 1 
k = sign(fi K ) = 0,±1 



i, = E(x l )cos<5 i cosa i 
y, = S(x») cos 5, sina t 
«4 = £(x0 sm< ^ 

r 7? cosh(x l /i?) fc = -1 
to< = < fc = 

{ R cos( X i/R) fc = +1 



k (X1X2 + yiV2 + ZjZ2) + W1W2 k = ±1 

a;ia;2 + 2/11/2 + 212:2 = 0. 



flcosh- 1 [(01,02) /R 2 ] k = -1 
X12 = ^ ^/ (a x - a 2 ,ai - a 2 ) fc = 

i?cos _1 [(ai,a 2 ) /i? 2 ] fc = +1. 



(33) 
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